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ABSTRACT 

We examine to what extent the inferred surface temperature of magnetars in quies- 
cence can constrain the presence of a superfluid in the neutron star core and the role 
of magnetic field decay in the core. By performing detailed simulations of neutron 
star cooling, we show that extremely strong heating from field decay in the core can- 
not produce the high observed surface temperatures nor delay the onset of neutron 
superfluidity in the core. We verify the results of Kaminker et al., namely that the 
high magnetar surface temperatures require heating in the neutron star crust, and 
crust heating is decoupled from cooling/heating in the core. Therefore, because crust 
heating masks core heating, it is not possible to conclude that magnetar cores are 
in a non-supcrfluid state purely from high surface temperatures. From our interior 
temperature evolutions and after accounting for proton superconductivity in the core, 
we find that neutron superfluidity in the core occurs less than a few hundred years 
after neutron star formation. This onset time is unaffected by heating due to core 
field decay at fields < 10^^ G. Thus all known neutron stars, including magnetars, 
without a core containing exotic particles, should have a core of superfluid neutrons 
and superconducting protons. 
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1 INTRODUCTION 

Neutron stars (NSs) begin their lives very hot (with tem- 
peratures T > 10^^ K) but cool rapidly through the emis- 
sion of neutrinos. Neutrino emission processes depend on 
uncertain physics at the supra-nuclear densi ties (p > Pnuc ~ 
2.8 X 10^^ g cm~^) of the NS core (se e iTsurutal 1 19981 : 
lYakovlev fc Pethick|[20o3 : |Page et al ]|2006l . for review) . Cur- 
rent theories indicate that the core may contain a neutron 
superfluid and proton superconductor or even exotic par- 
ticles, such as hyperons and deconfined quar ks (see, e.g., 
iLattimer fc PrakaslJ 12004 iHaensel et al.l |2007|. for review) 
The recent observation of rapid cooling ( Heinke fc Hdl2O10l : 



IShternin et al.|[201ll ) of the NS in the Cassiopeia A super- 
nova remnant provides the first constraints on the critical 
temperatures for the onset of superfluidity of core neutrons 
Tent (in the triplet state) and protons Tcp (in the singlet 
state), i.e., Tent ~ (5 - 9) x 10^ K and Tcp ~ (2 - 3) x 10^ K 
(|Page et al.ll201ll : IShternin et al.ll201ll ). 

Anomalous X-ray pulsars and soft gamma-ray repeaters 
form the magnetar class of NSs, i.e., systems which pos- 



sess superstrong magnetic flelds (B > 10^'' cj^, and their 
strong fie lds likely power the activ i ty seen in these ob- 
jects (see IWoods fc ThompsonI l2006l : iMereghettil |2008| . for 
review). One notable property of magnetars is that their 
observed surface temperatures Ts in quiescence are signif- 
icantly higher than those of other NSs of a similar age. 
In fact, they are too high for NSs that cool passively, 
i.e., without an additional source of internal heat (accre- 
tion heating can be excluded by, e.g., non-detections of 
binary companion or disk emission). An interesting prob- 
lem concerns the heat generated from magnetic field de- 
cay, which has been proposed to be the source for the high 
temperatures of magn e tars dThompson fc Duncan 19961: 
Hevl fc Kulkarnil 1 19981 : IColpi et all l2000l : lAguilera et all 



2008bl ). This heat can strongly influence the time/age 



at which the core becomes superfluid if heating/field 

decay occu rs in the core (iThqmpson fc Duncaiil Il996l : 

lArras et al.] |2004| : iDaU'Osso et al.l |2009| ). The problem is 
important since the presence of superfluid components 
has a strong impact on magnetar interior dynamics, such 
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^ However , there exists an apparently low surface magnetic field 
magnetar jRea et al.ll2010l) . 
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as th e mechanism for producing ghtches (see, e.R., Sauls 
198E ). fluid oscillations (see, e.g.. IPassamonti fc Andersson 



201 j). magnetic fleld and rotati onal evolution (see, e.g. 



Glampedakis fc AnderssonI 201l|). and magnetohydrody- 



produce surface temperatures that are high enough to ex- 
plain magncta r obse rva tions . 

Osso et al.l (|2009i) criticize the work of 



nami cal equilibrium ( Glampedakis et al]|2013 : lLander et al.l 
120121 ). 

We address this problem by conducting detailed calcula- 
tions of the thermal evolution of a NS with various prescrip- 
tions for an internal heat source that can be associated with 
magnetic field decay. We show that, regardless of the mag- 
netic field strength and detailed mechanism for field decay 
in the core, the heat generated by the decay is insufficient 
to power the surface emission of magnetars. Furthermore, 
by accounting for the effects of proton superconductivity 
on field decay, we find that this core heating is not strong 
enough to balance neutrino cooling and cannot delay the 
onset of core neutron superfluidity for flelds < lO'^'' G (on- 
set at age < a few x 100 yr). Thus the cores of all currently 
known magnetars should be in a superfluid state. 

We briefly describe past works and note key findings 
and assumptions made in these works that we improve upon 
here. lArras et al.l (|2004l ') consider coupled magnetic field de- 
cay and thermal evolution of magnetars. The internal field 
and superfiuid temperatures are assumed to be > 10^^ G 
and Tcp = 5 xlO^K and Tent = (5-9) x 10* K, respectively. 
lArras et all (|2004l ) find that magnetic field decay can delay 
the transition to core neutron superfiuidity and maintain a 
relatively high surface temperature to ages ~ 10'^ — 10^ yr, 
depending on Tent. However, their calculation only consid- 
ers volume-averaged quantities and thus assumes that the 
NS interior is isothermal (see Section [2. 2p . which cannot be 
the case if there is a localized heating source such as field 
decay in the crust or core due to the large but finite thermal 
conductivity of NS matter. 

iKaminker et al.l (|2006l ) calculate the evolution of the 
temperature profile [i.e., T{p)\ by solving the energy balance 
and flux equations [see eqs. ([2]) and ((3])]. They demonstrate 
that, in order to explain the observed high surface temper- 
atures, magnetars require a heat source and, most impor- 
tantl y, this heat source (e.g., from field decay; IPons et al.l 
l2009l ) must be located in the outer crust; if the heat source 
is located too deep in the NS interior (e.g., in the core), then 
neutrino emission efficiently removes the heat locally, and 
the surface temperature cannot be increased sufficiently to 
match the observed values. Furthermore, the outer crust is 
thermally decoupled from the core so that heating of the 
crust does not affect (neutrino) cooling of the core. Be- 
cause o f this thermal decoup ling, the results and conclu- 
sions of lKaminker et al] (|2006l ) are not particularly sensitive 
to the state of matter in the core, e.g., superfiuidity of the 
core nucleons. Though no quantitative results are shown, 
IKaminker et al.l (|2006l ) state that cooling calculations with 
the effects of inner crust superfluidity and neutrino emis- 
sion by Cooper pair formation produce different temperature 
proflles, but these effects d o not change the surfac e temper- 
ature. In follow-up work, IKaminker et al.l (|2009l ) improve 
their calculations and examine the effects of light-element 
accreted envelopes and anisotropic heat conduction due to 
the magnetic field in the envelope and outer crust. While 
accreted envelopes can give higher surface temperatures for 
the same core temperatures, core heating is still unable to 



IKaminker et al. I (|2006l ). arguing that the phenomeno- 
logical heating model considered is independent of magnetic 
field and decays on the wrong timescale [see eq. {TS))] 
and that the correct heating should depend inversely 
on temperature. This last point implies that more heat 
is generated at lower temperatures and can lead to an 
equilibrium between n eutrino cooling and h eating from 
field decay (see also iThompson fc DuncanI [T996. ) . The 
equilibrium temperature (a bove 10^ K) can be ma intained 
for ^ 10**yr. The resuhs of iDaU'Osso et al.1 (|2009l ') suggest 
that the core of magnetars are not superfluid until after 
this time, since the critical temperature for the onset of 
neutron superfluidity is (5 - 9) x 10* K (|Page et al.ll201ll : 
IShternin et al. 2011 ). However no cooling calculation is 
performed bv lDaU'Qsso et al.1 (120091 ). 

Here we perform NS cooling simulations to determine 
the role of core heating by magnetic field decay. We do not 
examine in detail field ev olution and heating in the crust 
(see, e.g.. |Pons et~al]|2009l). Rathe r we u s e the phenomeno- 
logical model of IKaminker et"all (|2006l . |2009| ) to demon- 
strate the effect of crust heating on the NS cooling behavior. 
This will be adequate for our purposes since the goal here 
is to assess the importance of core heating in magnetars. 
As we will show, as long as the magnetic field decay time 
is longer than the cooling time, the core heating we use is 
the maximum one can use. While this maximal heating can 
delay the on set of neutron superfiuidity [o n the timescale 
estimated by IThompson fc DuncanI (|l996l ) ; iDall'Osso et al.l 
(|2009l )]. it cannot reproduce the high surface temperatures 
seen for magnetars. We arrive at this same conclusion us- 
ing the prescription f or heating and field decay given in 
iDaU'Qsso et"alj |200i). Furthermore, the above works ne- 
glect the effects of proton superconductivity; when this is 
taken into account, there is no delay in superfluidity on- 
set for any reasonable core magnetic field strength. In Sec- 
tion (21 we describe the thermal evolution equations and in- 
put physics, including superfluid properties and prescrip- 
tions for internal heat sources. Section [3] presents the results 
of our calculations. We summarize our results and discuss 
their implications in Section 31 



2 NEUTRON STAR COOLING MODEL 

2.1 Core and Crust Composition 

We use a stellar model based on the Akmal-Pa ndharipande- 
Ravenhall (APR) equation of state (E OS) (jAkmal et al.l 
Il998l: iHeiselberg fc Hiorth- Jensen! [l999l ). specifically APR 
fl (|Gusakov et al.ll2005l ). For a NS mass M = 1.4 M©, the 
NS radius is R = 12.1 km, and the central density is 



Pc 



9.47 X lO^'' g cm ^. The crust compositioi|f| (in par- 



2 The maximum NS mass for this EOS is 1.9 23 Mp,, which is 
below the highest measured mass of 1.97 l lDemorest et al] 
I2OIOI) . Nevertheless, APR is typical of EOSs that yield higher 
maximum masses, and our conclusions do not depend strongly 
on the specific EOS. 

^ See [http: / / relativity.livingreviews.org/ Articles/lrr-2008-10/| 
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ticular, charge number Z and mass number A of the ions) 
as a function of density p is deter mined by interpola ting 
the values taken from Table VII of iRiister et all (j2006l ) for 
8. 02 X 10^ K cm~^ - P.- ^-. ^'^ ^ g cn 
of lNegele fc Vautherini (| 19731 ') for 4.67 x 10 
7.94 X 10^ ^ g cm ~'^ and estimated from Figs. 6 and 9 of 
lOvamatsul (| 19931 ) for p = (1, 1.5) x 10^* g cm~^ We note 
that the melting ternperat ure Tmeit is given by (see, e.g., 
IShapiro fc Teukolskvl[l983 ) 



^ and Table 3 

g cm^^ < /5 < 



3.04 X 10^ K (Z/26)''''^(170/rn0a;, 



(1) 



where x = [{Z/A){p/W^ g cm"^ 



and Tm ~ 170 is the 
value of the Coulomb parameter Fc [= Z^e^/akT, where 
a = (3/47rni)^'''' and ni is the ion number density] at which 
melting occurs. 



2.2 Equations for Neutron Star Cooling 

The evolution of the interior temperature T{r,t) of an iso- 
lated NS is determined by the relativistic equations of energy 
balance and heat flux, i.e.. 



A-Kr^ dr \ '^J dt " 



4nr^ d\ 



(2) 
(3) 



where Lr is the luminosity at radius r, C is the heat ca- 
pacity, is the neutrino emissivity, Qh is the internal 
heati ng source, and K is the thermal condu ctivity (see, 
e.g., iThornd [l977l : IShapiro fc Teukolskvl [19831 '). The met- 
ric function A is defined by the enclosed mass m{r), i.e.. 



= [l — 2Gm{r)/c^r\ ^^'^ , $(r) is the metric function 



eA(.) 

(gravitational potential in Newtonian limit) that is obtained 
from d^/dr = -{dP / dr) /[p{l + P/p)], and_e2^1= e*'^^ 
at the NS surface. We define (see, e.g.. Ivan Ripei|[l99ll ') 

f = e'^T and Z = e^*L„ (4) 

so that eqs. ^ and ([HJ become 



or at 



. 2 A+2* I 

- Anr e (Sv 



Oh) 



L = -4^r27^e-*+* 



dT 
dr ' 



(5) 



(6) 



respectively. Thus the inputs that determine how a NS 
cools are the heat capacity C, neutrino emissivity e^, ther- 
mal conductivity K, and heat source Qh- The following 
sections discuss the physics that determine the four phys- 
ical parameters, C, e^, K, and Qh- We note that, due 
to the high thermal conductivity, NSs without additional 
heat sources (Qh = 0) become essentially isothermal af- 
ter ^ 10 - 100 yr (iLattimer et al.lll994l : ICnedin et al]|200ll : 
lYakovlev et al.ll201ll ). When isothermal, the thermal evolu- 
tion is determined solely by eq. ((Sj, which can be simplified 
by integrating over the volume. The evolution of the interior 
te mperature is s o lved b y a metho d similar to that described 
in iGnedin et all (|200ll ) (see also lYakovlev fc Pethick|[200i : 
I Page et al.l 200t , for review) . 



Table 1 . Phcnomcnological parameters for superfluid gap energy 









shallow 


deep 




proton 


neutron 


neutron 


neutron 




singlet 


singlet 


triplet 


triplet 


Ao (MeV) 


120 


68 


0.068 


0.15 


ko (fm-i) 





0.1 


1.28 


2 


fci (fm-2) 


9 


4 


0.1 


0.1 


k2 (fm-i) 


1.3 


1.7 


2.37 


3.1 


k-i (fm-2) 


1.8 


4 


0.02 


0.02 



2.3 Superfluid Properties 

Superfluidity has two important effects on neutrino emission 
and NS cooling: (1) suppression of heat capacities and emis- 
sion mechanisms, like modifled Urea processes, that involve 
superfluid constituents and (2) enhanced emission due to 
Cooper pairing of nucleons when the temperature decreases 
just below the c ritical value (see lYakovlev fc Pethickll2004l : 
IPage et al.ll2006l . for review). 

The critical temperatures for superfluidity are approx- 
imately related to the superfluid energy gap A by kTc « 
0.5669A for the singlet (isotropic pairing) gap and kTc ~ 
0.8418A for the maximum triplet (anisotropic pairing) gap. 
We use the parametrization for the gap energy 



A(fcFx) = Ao 



(fcpx - fco)^ 



(fcp 



(fcpx - ko)^ + ki (fcpx - ^2)2 -I- fcs ' 



(7) 



where ppx — ftfcpx — h{3Tv^n,c)^^''^ and rix are the Fermi 
momentum and number density, respectively, for particle x 
and Ao, ko, ki, fe, and fca are flt parameters for particu - 
lar supe rfluid gap models given in [Andersson et al.l (l2005h 
(see also iLombardo fc Schu£3l200ll : [Kaminker et al]|2002l '). 
For singlet neutrons in t he crust, we use model a, which 
represents the results of IWambach efsH l)l993h . For sin- 
glet proton s in th e core, we use a model tha t is similar to 
IChen et al.l (|l993l ) (see also lPage et"aLll2004h . We consider 
two cases for triplet neutron pairing in the core. Either the 
pa iring is the "shallow" model /, which represents the results 
of lElgar0v et al.1 (|l996l ). or the pairing is a "deep" model. 
The deep model is a A that produces a rcnt(p) similar to 
what is needed to fit the observed surface temperature of the 
neutron s tar in Cassiopeia A and neutron stars that are ol d 
and hot (|Gusakov et al.1 12004| . l2005l : IShternin et al.1 [2OI1I '). 
The parameters for the superfiuid gap energies are given in 
Table [1] and the gap models are illustrated in Figs. [1] and (2] 



2.4 Heat Capacity 

In the core, the total heat capacity is the sum of partial heat 
capacities due to neutrons, protons, electrons, and muons. 
In the crust, the total heat capacity is the sum of partial 
heat capacities due to free neutrons, ions, and electrons. 
The partial heat capacity for a strongly degenerate fermion 
particle species x is 



Cx 



m*ppxfc T 
3fi3 ' 



(8) 



where m* is the effective mass. We assume mj = 0.7 for 
neutrons and protons and take = {m^c^ + pp^c^y^^ for 
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Figure 1. Supcrfluid gap energies (left: singlet; right: triplet) 
as a function of Fermi wavenumber for neutrons fcp„ and pro- 
tons fcpp. The neutron singlet (dot-short-dashed) model is from 
IWambach et al.l lll993l'). the proton singlet (dashed) model is sim- 
ilar to TcheiT^t al. I l ll993l) . and the neutron tri plet (dot-long- 
dashe d) model is either a shall ow model from [Elgargy et alj 
lll996l ') or a deep model similar to IShternin et al.l boill) . 



electrons and muons. A reduction of the heat capacity due 
to nucleon superfluidity can be taken into account by using 

SF 



a 



(9) 



where the reduction factor i?^ only depends on the ratio 
T /Tex and whe ther the supcrfluid is o f singlet or triple t type 
l|Levenfish fc Y akovlcv 1994a: Yako vlev et al.l llQQQbl ). The 
ion heat capacity is (jvan 



(3/2)A: X <^ 2/d(T/Td) 



log 150 



if Tc < 1 

if 1< Tc < 150 

if Tc > 150 



(10) 



where 

r o.87r*(r/rD)' 

foiT/Tu) = <j 1 - 0.05(Td/T)2 

[ 1.70(T/Td) + 0.0083 

and the Debye temperature is 

Td = 3.48 X 10^ K (ZM)(p/10^ g cm-3)i/2. 

2.5 Neutrino Emissivity 



if T/Td < 0.15 
if T/Td > 4 (11) 
otherwise 



(12) 



For the NS core, we calculate neutrino emission due to the 
neutron and proton branches of the modified Urea process 
and neutron-neutron, neutron-proton, and proton-proton 
bremss tr ahlung using em issivities from I Yakovlev et al.l 
l|l999bl ): |Page et al.1 (|2004l ). When neutrons and/or protons 
are su perfluid, we account for suppression of the above pro - 
cesses (jLevenfish fc Yakovlevlll994bl : lYakovlev et~al]|l999bl ) 
and neutrino emi ssion due to Cooper pairing of the sup er- 
fluid component ([Yakovlev et al.lll999bl : [Page et al.ll2009l ). 



10'» 
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core 
T 
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Figure 2. Superfluid critical temperatures as a function of den- 
sity (top panel) and normalized stellar radius (bottom panel). 
Critical temperatures Tens (dot-short-dashed). Tent (dot-long- 
dashed), and Tcp (dashed) are for neutron singlet, neutron triplet, 
and proton singlet, respectively. Neutron triplet pairing is taken 
to be described by either the shallow or deep model (see text). 
Vertical solid lines indicate the approximate boundaries be- 
tween the core and inner crust (at near nuclear saturation, i.e., 
P ~ Pnuc/2) and inner and outer crusts (at neutron drip, i.e., 
p 4 X 10'^'^ g cm~^). 



For the APR I EOS (see Section I2.1|l , neutrino emis- 
sion by the direct Urea process occurs when M > 1.829 
(cent ral densities above 1.680 x 10^^ g cm~'^: [Gusakov et al.l 
l2005h : for this present work in which we consider only M = 
1.4 Mq, we will negle ct this process (as well as other 'fast' 
cooling processes; see lYakovlev fc Pethicklliooil ; IPage et al.l 



l2006h . As shown by iKaminker et al. ( 20061 ). cooling by dt 



rect Urea has no effect on the surface temperature because 
of thermal decoupling between the core and outer crust. In 
addition, since direct Urea is more effective at cooling the 
core, it would speed up the onset of superfluidity. Hence its 
inclusion would make our conclusions even more robust. 

In the crust, we account for electron- nucleon, neutron- 
neutron, and neutron-nucleon bremsstrahlung, plasmon de- 
cay, and e'-e" *" pair annihila tion using emissivities from 
lYakovlev et al.1 (|l999al . l200lh When neutrons are super- 
fluid (in the singlet state), we account for the suppression 
of neutron-neutron and ne utron-nucleon bremsstrahlung 
(lYakovlev et al.lll999bl . 1200 j ), as well as inc luding neutrino 
emissi o n due to neutro n Cooper pairing ([Yakovlev et al.l 
Il999bl ; [Page et al.|[2009l ). We neglect neutrino synchrotron 
emission in a magnetic fleld since this additional emis- 
sion would cause ev en more rapid cooling of the crust 
(lYakovlev et al.ll2001^ ■ 



2.6 Thermal Conductivity 

For the core, we sum t hermal condu ctivities d ue to neutrons, 
electrons, and muons (jFlowers fc Itoh , 19791 , [l98ll ). We use 
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the results of iBaiko et alJ (|200 j ) to calculate the neutron 
thermal conductivity, accounting for neutron-neutron and 
neutron-proto n collisions and superfluid su ppression. We use 
the results of IShternin fc YakovlevI (|2007l ) to calculate the 
electron and muon thermal conductivities and account for 
superfluid suppression. For the crust, we use CONDUCTofl 
which implements the latest advancements i n calculating 
thermal conductivities JPotekhin et al.l 1 19991 : ICassisi et al.l 
I2OO7I : IChugunov fc Haenselll2007l ). We assume no contribu- 
tion by impurity scattering, which only becomes important 
at high densities and low temperatures. 

At high magnetic fields and low temperatures, electron 
mot ion is strongly influenced b y the magnetic field (see, 
e.g.. lYakovlev fc Kaminkeilll994l ). As a result, the electrical 
and thermal conductivities become anisotropic, depending 
on whether electro ns move paralle l or transverse to the di- 
rection of the field (|Potekhinlll999l '). However, electron mag- 
netization and anisotrop ic heat co nduction are weaker at 
higher temperatures (see lLaill200ll . for review). In particu- 
lar, magnetic field effects are minimal when the field is in 
the non- or weakly-quantizing regime, i.e., when T > Tb 
and p ^ /5b, where 



(13) 
7.046 X 



Tb = 1.34 X 10'" K (B/IO'* G)(l + a;^)"'/^ 

X is defined below eq. ([T]), and pB 
10^ g cm-^(A/Z)(_B/10^'' G)^/^ Since we are primarily 
concerned with effects in the high density core of young 
NSs with high interior temperatures and t he core is effec- 
tively therm ally decoupled from the crust (|Kaminker et al.l 
I2OO6I . |2009| ), we neglect magnetic field effects on the con- 



ductivities below the heat blanket (see also Section 12. 7|) . 
In fact, recent works show t hat the effects of superfluid 
phonons ("Aguilcr a et al.ll2009l ') and toroidal magnetic fields 
IjPons ct al. 2009.') in the crust can very effectively smooth 
out temperature variations induced by anisotropic heat con- 
duction. 



2.7 Envelope 

The outer layers (envelope) of the NS crust serve as a 
heat blanket, and there can exist a large temperature gra- 
dient between the bottom of the envelope (at ponv and 
Tenv ) and the surface. \Ve assume pen v — 10" 



g cm ^ (see , 

e..g, iGudm undsson ct al.l ll982l : lYakovlev fc Pcthickl [20o3 ). 
IChang fc Bildstcn (200i ) (see also Chang ct al. 20jJ) show 
that, at the high temperatures present in magnetars, nu- 
clear burning very rapidly removes any surface hydrogen; 
nevertheless, for generality, we consider surfaces composed 
of either iron or light-elements. For iron, the relationship be- 
tween the temperature at the bottom of the heat-blanketing 
envelope Ten v and the effective tempe rature of the photo- 
sphere Ts is (|Gudmundsson et al.lll982^ 



Ts = 8.7 X lO"" K gl^,\ (Tenv/lO* K) 



11/20 



(14) 



where (7s,i4 

j2 D^-l/2 



Ps/lO cm s and gs = (1 ^ 
2GM/c^R)-'''-'GM/R^ is the surface gravity (ps,i4 = 1.55 
for our assumed NS model). Non-magnetized envelopes com- 
posed of light elements (H or He) have been considered in 



http: / /www. ioffo.ru/astro/conduct/ 



IPotekhin et~al] |T993) [Kaminker et al.l (|2009l) find a neg- 
ligible difference in t he results using eith er H or He at 
cnv > 10*^ K; see also lYakovlev et~aLll201ll . for carbon en- 
velopes]. Due to the higher thermal conductivity (oc Z^^) 
of light elements, the difference between Ts and Tenv is 
larger for a light-element envelope compared to an iron en- 
velope, i.e., Ts is higher for a given T^nv and Tenv is lower 
for a given Ts. For a fully-accreted light element envelope 
(|Potekhin et al.lll997l '). 



1.43 X lO'' K gl'-^^ (Tenv/10« K) 



17/28 



(15) 



For clarity, we show the above relations [eqs. (|14|) and (|15|l ]: 
however, we u s e the more accurate expressions given in 
IPotekhin et~al] (120031 ) in our calculations. 

Magnetized envelopes produce higher (lower) Ts for 
a given Tenv due to enhanced (reduced) thermal conduc- 
tion along (across) the magne tic field, and the effect is 
stronger for accreted envelopes ([Po tckhin fc Ya kovlevll200ll : 
IPotekhin et~al]|2003t ). iKaminkeret al. (2006, 2 0"09l 'l include 
the effects of magnetic fields (^ 10^'' — lO^'' G) in the enve- 
lope and inner and outer crusts and find that, though the 
temperature profile in the crust is modified, the average sur- 
face temperature is only weakly affected by anisotropic heat 
conduction. We mostly neglect magnetized envelopes since 
their effects would not significantly affect our conclusions 
(see Section [3. 5|) . 



2.8 Heat Source 

We consider several prescriptions for an internal heat source. 
For heating in the cr ust, we use a model similar to the phe- 
nomenological one of iKaminker et al] (|2006l . [2"009t ). i.e.. 



Qh = Qo exp{-[(p - ph)/Ah]-'}exp(-t/rh); 



(16) 



the exact form of heatin g is u nimportant, as demon- 
strated bv iKaminker et al.l (|2009h . We take Th = lO"* yr, 
Qo = 3 X 10^° ergs cm-3 s"\ ph = 6 x 10^° g cm-^ 
and Ah = 2 X 10"° g cm"''. The total heat luminosity 
Th = [Qhe^* Anr'^dr X lO^J ergs s"" at t < Th. Note that 
IKaminker et al.l (|2006l . |2009| ) use a heating decay timescale 
Th = 5 X lO** yr; the exact value is unimportant for our pur- 
poses, as long as it is approximately the age of the oldest 
(high surface temperature) magnetar, since Th is approxi- 
mately the time during which a high surface temperature 
can be maintained by this crust heating (see Section [3. 5p . 

In the core, heating is supplied by magn etic field decay 
due to ambipolar diffusion (jGoldreich fc ReiseneggGriil992i '). 
We estimate a heating rate 



B dB 



pnuc 



B^ 



-2/3 



1 X 10^° ergs cm"^ s"" 



V109 K 



B 



1016 G 



(17) 



where the ambipolar diffusion timescale r^'"'' is given by 
(|Goldreich fc Reiseneggeilll992l ) 



2.5 X 10' yr (T/1 km)"(p/pnnc)'''' 
x(T/10^ K)''(B/10''^ G)"' 



(18) 



and we neglect here the proton fraction dependence and as- 
sume the field decay lengthscale T (= 1 km) is constant for 
simplicity. Heating is taken to occur throughout the core and 
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is thus very st ronR. Note that eg. (1171) is identical to what is 
considered bv lPaH'Osso et ahl (|2009l ). except for a small dif- 
ference in the numerical coefficient. It is also worth noting 
that r^™*^ is the timescale of the "solenoidal" mode. This 
mode is controlled by collisions between the different parti- 
cle species in NS matter and operates at relatively high tem- 
peratures (> a few X 10* K). Superfluidity can significantly 
modify r™*^ and thus the heating rate given by eq. (|17|) . 
This is discussed in more detail in Section [3.61 

It is not well-understood how magnetic field decays (if 



at all) in the core at times t < 10 yr ( see. e.g..lHaensel et al 



1990'; Goldrcich fc Rciscncgi^e3 ll992l : [Thompson fc Duncan 
1996. : .Glampedakis ct al.. ,201lh . Therefore, to describe the 



decay of the core magnetic field, we use the simple formula 
Bo 



B{t) 



1 + t/rac 



(19) 



where rdccay is the field decay timescale, which we choose to 
be 10^ yr. Though eq. (|19|) has no physical basis, it has a 
form similar to ones considered in the literature, usually with 
the denominator (i.e., time-dependence) to some low-order 



power (see, e.g., Thompson fc DuncanI 1 19961 : IColpi et"al] 



I2OQOI : iDaU'Osso et al.ll2009l )l^ However. recaU that 



our aim 

is to obtain a maximum core heating rate, in order to test 
the extent to which core heating can produce the observed 
(high) surface temperature of magnetars and delay the on- 
set of superfiuidity. It is clear that eq. fTTJ is maximal be- 
fore significant field decay, i.e., B = Bo at t <^ Tdccay. [At 
t = 10^ yr, we find Lh ~ 2 x 10^* ergs s-^(B/10^^ G)''.] 
Therefore, as long as Tdecay is longer than the thermal dif- 
fusion timescale, this core heating will have a maximum 
effect on the surface temperature. We also conduct simu- 
latio ns utilizing the heatin g rate and field decay formulae 
from lDall'Osso et all (|2009h and find very similar results at 
the ages of interest here (< 10^ yr). 



3 RESULTS 

We assume a constant Te* = 10^" K at t = 0. For clarity, 
the superfluid results shown in all subsequent figures use the 
shallow neutron triplet model (see Section [2.3p : the super- 
fiuid results using the deep neutron triplet model are quali- 
tatively similar in the age range considered here (< 10^ yr). 
Also we assume an iron envelope, unless noted otherwise. 
For models which include core heating, we primarily assume 
an initial magnetic field Bo = 10^^ G, though we also con- 
sider lower, more realistic initial field strengths, as well as 
extreme fields. 



^ iTIiompson fc DuncanI l l 19961 ) lias a 1/14-th power of 
the denomin ator and a de c ay ti mescale proportional to 
.^amb^ while iDaU'Osso et al.l ^20091') has a 5/6-th power 
of the denominator and a decay timescale given by 
1.03 X 10* yr (B/IO^^ G)-^/5(p/10i5 g cni-3)2/5 r^^ie dif- 
ferent evolution laws reflect the different ambipolar diffusion 
modes (solenoidal vs non-solenoidal) considered by these authors. 
Note that these decay timescales are shorter for higher magnetic 
fields. Incidentally, numerical simulations of field evolution in 
the crust yield eq. I I19I I. with Tdecay given by the timescale for 
Hall drift at times mu ch shorter than the Ohmic decay timescale 
jAguilera et al.ll2008a| [bll. 



10" 



'10" 



CD 



108 



outer crust 


. 1 . 11 1 1 1 1 1 . 1 

inner crust 
: 10-4 yr^^^^-^.^J 


core 

- 










" — ^-^-..^^ 






100 h- 












10= 












lO' yr 












10 



11 



12 13 
log p (g 



14 



15 



Figure 3. Temperature profiles for cooling model with no super- 
fiuidity and no heating (dotted lines). The six profiles are at ages 
t = 10"* (top), 0.01, 1, 100, 10^, and 10* yr (bottom). Thin solid 
lines denote the melting temperature Tmelt, Debye temperature 
Td, and Tb, where Tb is calculated assuming B = 10^^ G [see 
eq. I|13|l ]. Vertical solid lines indicate boundaries between core and 
inner crust and inner and outer crusts. 



3.1 Cooling with no superfluid nor heating 

Figure [3] shows the temperature as a function of density at 
different ages for a NS cooling model that neglects super- 
fiuidity and has no additional sources of internal heating 
(Qh = 0). The se temperatur e profi les are very similar to 
those shown in iGnedin et al.l (|200ll ). At very early times, 
the core cools more rapidly than the crust via stronger neu- 
trino emission, so that the crust is generally at higher tem- 
peratures. A cooling wave travels from the core to the sur- 
face, bringing the NS to a relaxed, isothermal state. De- 
pending on the p roperties of the crust , the relaxation time 
is ^ 10 - 100 yr (iLattimer et al.l 1 19941 : iGnedin et al] I2OOII : 
lYakovlev e^aL[|2011^ . Formation of the inner and outer 
crusts begins at ~ 1 hr and ~ 1 day, respectively, and is 
mostly c omplete after ~ 1 mon th and ~ 1 yr, respectively 
(see also lAguilera et al.l l2008al ) . The temperature profiles 
demonstrate the need to account for magnetic field effects 
in more accurate models of the crust and envelope when 
T <Tb and temperature gradients are significant (see Sec- 
tions [M] and [gjl) . 



3.2 Cooling with np superfluid and no heating 

Figure |4] shows T{p) at different ages for a cooling model 
that neglects heating but includes superfluidity of neutrons 
and protons, with critical temperatures also shown. When 
superfiuidity is taken into account, we see the impact of 
the two effects mentioned in Section 12.31 Slower cooling in 
the core after protons become superconducting and faster 
cooling after neutrons become superfluid due to neutrino 
emission from Cooper pair formation. The latter is strongest 
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Figure 4. Temperature profiles for cooling model with superflu- 
idity and no heating (long-dashed lines). The six proflles are at 
ages t = IQ-* (top), 1, 100, 500, 10^, and 10'' yr (bottom). Also 
plotted for comparison are profiles (dotted lines) for the model 
with no superfluidity at ages t = 10~*, 1, and 100 yr (see Fig.[3ll. 
Critical temperatures for neutron singlet (dot-short-dashed) , neu- 
tron triplet (dot-long-dashed), and proton singlet (short-long- 
dashed) are shown. Vertical solid lines indicate boundaries be- 
tween core and inner crust and inner and outer crusts. 



in regions near the critical temperature. Proton supercon- 
ductivity occurs at ~ 1 min, and much of the core is su- 
perconducting after ~ 1 yr. Core neutrons start becoming 
superfluid at around a few x 100 yr. 

3.3 Cooling with superfluid and crust heating 

Figure [5] shows temperature profiles for a cooling model 
that includes superfluidity and crust heating [see eq. (|16[l ]. 
T he proflles with c r ust h e ating are similar to those shown 
in iKaminker et al.l (|2006l . |2009| ). In particular, it is clear 
that an additional heat source in the outer crust can very 
effectively maintain a high temperature near the NS sur- 
face [with redshifted surface temperature > 3 x 10® K, 
where T;r° = rs(l - 2GM/c^RY^'^, for times longer than 
Th (= 10^ yr)], but this strong heating in the crust does 
not prevent the core from cooling rapidly. Due to thermal 
decoupling between the outer crust and core, the core tem- 
perature drops below the critical temperature for neutron 
triplet superfluidity at a few x 100 yr, whether or not there 
is crust heating. 

3.4 Cooling with superfluid and normal core 
heating 

Figure |6] shows temperature proflles for a model that in- 
cludes superfluidity and core heating due to magnetic fleld 
decay with an initial field Bo = 10^® G and Tdccay ~ 10* yr 
[see eqs. (|17p and (|19|) ]. We see that, with such strong heat- 
ing, the core temperature stays above the critical temper- 



Figure 5. Temperature proflles for cooling model with super- 
fluidity, crust heating, and heating timescale rh = 10* yr (short- 
dashed lines). The six profiles are at ages t = 10~* (top), 1, 100, 
500, 10^, and 10* yr (bottom). Critical temperatures for neutron 
singlet (dot-short-dashed) , neutron triplet (dot-long-dashed) , and 
proton singlet (short-long-dashed) are shown. Vertical solid lines 
indicate boundaries between core and inner crust and inner and 
outer crusts. 



ature Tent for t 2> 100 yr. However, the extra heat created 
in the core from field decay is efficiently removed by core 
neutrino emiss ion processes, with their strong temperatur e 
dependencies (jYakovlev fc Pethickl [iooi : IPage et al.ll2006l ): 
while the temperatures in the crust are higher than those 
from models without any heating, they are significantly 
lower than those from models with crust heating. 

3.5 Comparison with magnetar surface 
temperatures 

Figure [7] shows the evolution of redshifted surface temper- 
atures for the cooling models plotted in Figs. |3] - |6l 
For models with crust or core heating, we also show cool- 
ing curves with fully accreted light element envelopes and 
a surface magnetic field 10^^ G that is directed along the 
radial direction, which produce higher surface temperatures 
for the same outer crust temperature (jPotekhin et al.ll2003l : 
see also Section I2.7|l . The cooling curves are compared to 
the estimated ages and measured surface temperatures of 
magnetarjl and other neutron stars (see [Chcvalicr 200 ^ ; 
Yakovlev et~all l2008l : IHo fc Heink3 l2009l : licaminker etldj 
20091 , and references therein). Our results clearly indicate 

S Data are taken from the McGill SGR/AXP Online Catalog 
at 'http: / /www. physics. mcgill.ca/~pulsar/magnetar/main. html 
Magnetar ages and surface temperatures are not well-determined 
due to additiona l systematic uncertainties that are large (see, e.g. , 
IKaminker et al.i r2006. 200i, for details). IKaminker et al.l ||2009| ') 
consider the more reliable surface luminosity. For simplicity, we 
consider the "magnetar box" illustrated in Fig. [7] 
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Figure 6. Temperature profiles for cooling model with superflu- 
idity, normal core heating, initial magnetic fleld Bq = 10^^ G, and 
field decay timescalc T^ccay = 10^ yr (solid lines). The seven pro- 
files are at ages i = lO^"* (top), 1, 10^, 10* yr, 3 X 10*, 5 X 10*, 
and 10^ yr (bottom). Critical temperatures for neutron singlet 
(dot-short-dashed), neutron triplet (dot-long-dashed), and pro- 
ton singlet (short-long-dashed) are shown. Vertical solid lines in- 
dicate boundaries between core and inner crust and inner and 
outer crusts. 



that heating in the (outer) crust can produce surface temper- 
atures which match the high magnetar surfa ce temperatures; 
this is in a greement with the findings of iKaminker et al] 
(|2006l . l2009l ). On the other hand, even using a more real- 
istic core heating prescription that includes magnetic field 
strength and temper ature dependencies (as suggested by 
iDairOsso et al.ll2009l ). heating of the core by magnetic field 
decay produces that are too low to explain the observed 
temperatures of magnetars. Only for extreme (_Bo > 10^® G) 
and long-lived (rdccay -SJ 10'' yr) magnetic fields and a mag- 
netic, fully accreted light element envelope (which is un- 
hkely due to diffusive nuclear burning in the hot surface 
layers; IChang fc Bildstenll2004l ) could a surface temperature 
at the low end of magnetar temperatures {T^ ~ 3 x 10^ K) 
be produced. Therefore core heating alone cannot explain 
magnetar surface temperatures. 

3.6 Onset of core neutron superfluidity with 
proton superconductivity 

From Figs. |4] - [6l we see that the onset of core neutron 
superfiuidity can be delayed from age a few x 100 yr when 
the initial core field Bo < a few x 10^* G to, e.g., ^ 10* yr 
when Bo = 5 X 10^^ G (and rdocay = 10* yr) or 3 x 10* yr 
when Bo — 8x 10^® G (and Tdecay ~ 10^ yr). This is based on 
our simple model of core heating, especially the assumption 
that field decay occurs on the timescale given by eq. (|18p . For 
the deep neutron triplet model (see Sec. 12. 3p . onset occurs 
earlier since the inner core cools faster and the maximum 
Tent is higher than in the shallow model (see Fig. [2]). 



magnetars 




t (yr) 



Figure 7. Rodshifted surface temperature evolution for mod- 
els with no superfluid (dotted line), superfluid and no heat- 
ing (long-dashed line), superfluid and crust heating for heating 
timescale = 10* yr (short-dashed lines), and superfluid and 
core heating for initial magnetic field Bq = 10^® G and field de- 
cay timescale Tdocay = 10^ yr (solid lines); the thin upper lines 
have fully accreted light-element envelopes (with a 10^^ G ra- 
dial surface magnetic field), while the thick lower lines have iron 
envelopes. Data points are magnetars and other neutron stars 
ta ken from the McGill SGR/AXP Online Catalog and those liste d 
inlChevalieil ll2q05l) : pYakovlev et al.l ||2008|) : IHo fc Heinkj | |2009^ : 
iKaminker et al .1 1|2009^ . respectively. 



The decay timescale given by eq. (|18|l is the timescale 
for the solenoidal mode of ambi polar diffusion of normal 
(non-superconducting) protons (|Goldreich fc Reiseneggeil 
1 19921 ; iThompson &i DuncanI [T996I ) . At lower temperatures 
(later times), there is a transition of ambipolar diffusion 
from the solenoidal to the non-solenoidal mode. This latter 
mode is driven by de parture from beta chemical equ ilibrium 
and has a timescale (jGoldreich fc Reiseneggerlll992l ) 

r"=°' ~ 240 yr {p/ p^^,f'\B G)-\T/l(f K)-\ (20) 

and the transition occurs at temperature 

T*'- = 7x 10* K (p/pn.c)'''''. (21) 

Since Tcp ~ (2 - 3) x lO'-* K l|Page et al.l 1201 ll ; 
IShternin et al.|[201ll ). proton superconductivity occurs very 
early. When protons are superconducting (and neutrons are 
normal), processes invol ving protons are suppressed (see 
iGlampedakis et al.ll201l] . for more detailed calculation). In 
particular, the beta equilibrium reaction rate is reduced 
by a factor RpA, w hich depends only on the ratio T/Tcp 
(jHaensel et al.ll200lh . and the solenoidal and non-solenoidal 
ambipolar diffusion timescales are given by 

rlf = {RpB/H,^)r^''* (22) 

r,r' = (B/i/ciJ?pA)r"=°', (23) 

respectively, where Hci {~ 10^^ G) is the first criti- 
cal magnetic field, Rp [— Rp + {m* /nie){Den/ Dp„)] is 

© 2012 RAS, MNRAS OGO.ITlfTTI 



Superfluid core in hot magnetars 9 



modified to account for electron- 
p is ttie r eduction fac to r of proton- 



and 



the reduction factor 
neutron collisions, Rp 
neutron collisions given 

I>en = 10^* 

4.8 X 10^* s^^ (p/pnuc)"^''^(r/10^ K)^ are the inverse re- 
laxation timescales for e-n a nd p-n collisions, respectively 
lYakovlev fc Shalvbkovlll990l ). The transition temperature 
also increases, i.e., 



m iBaiko et al] (j200lh . 

s ^ (/9/pn,e)-''"(r/10« K)2 and Dp„ 



i^c = (i?pAi?p)-'/'r"- 



(24) 



Thus the heating rate in the core due to magnetic field decay 
by ambipolar diffusion of normal or superfluid protons is 
given by eq. inillll where 



if T > Tcp 

if T^^ <T <Tc, 

if T < T''- 



(25) 



and illustrated in Fig. [S] for typical parameter values. Note 
the removal of one power of B in T^^^ at T < Tcp and the 
strong increase in t^^^ (oc R^^T~^) at low temperatures. 
Using this heating rate (that accounts for superconducting 
protons) in our cooling calculations, we find that there is no 
delay in the onset of core neutron superfluidity. Onset oc- 
curs at a fewx 100 yr, and the temperature proflle/evolution 
is very similar to the case without core heating (see Figs. [4] 
and[S| for any fleld < 10^® G (however, at the hi gh fleld end, 
super conductivity is expected to be suppressed: iBavm et al.l 
1 19691 ). This is because proton superconductivity greatly in- 
creases the timescale for fleld decay, which consequently re- 
duces Qh and renders core heating ineffective against cooling 
by strong neutrino emission. 
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Figure 8. Core heating and field decay timescales used in this 
work as a function of temperature for typical parameter val- 
ues. Long-dashed line is (solenoidal) ambipolar diffusion timescale 
^amh normal matter [see eq. I)18|l]. dot-dashed and solid lines 
are ambipolar diffusion timescales rJJ,™'' which account for proton 
superconductivity at Tcp = 3 X 10^ K and 6 X 10® K, respectively 
[see eq. 12511 ]. and short-dashed line is the field decay timescale 
"""decay [= 10* Y^'t ^^e cq. ||19| |] . Vertical dotted lines denote the 
critical temperatures for proton superconductivity Tcp and tem- 
peratures for transition between solenoidal and non-solenoidal 
modes of ambipolar diffusion T^^ [see eq. II24| |]. 



4 DISCUSSION 

We performed detailed calculations of neutron star cooling, 
including the effects of superfluidity and additional heating 
(due to magnetic fleld decay) in the crust and core. We flnd 
that magnetic field decay in the neutron star core cannot 
be the sole source powering the high observed surface tem- 
perature of magnetars in quie scence; the high temperature s 
require outer crust heating (jKaminker et al.l |2006| . [20091 ). 
Because of crust heating and effective thermal decoupling 
between the outer crust and core, the state of matter in 
the core cannot be deduced from these surface temperature 
measurements. By computing the evolution of the tempera- 
ture profile T{p), we determine the time when core neutrons 
first become superfluid, i.e., when T < Tcnt(p). We find that 
heating by field decay in the core (with fields < 10^^ G) can- 
not balance neutrino cooling and thu s cannot maintain rel- 
ative l y high core temperat ures (c.f., [Thompson fc DuncanI 
ll996l : iDaU'Osso et al.lf2009l '). As a result, onset of superfiu- 
idity for neutrons in the core cannot be delayed, and neu- 
tron stars possess superfiuid and superconducting cores af- 
ter a few hundred years; this does not strongly depend on 
the nuclear EOS. Since core heating is not significant, the 
temperature profiles T(p) and surface temperature evolution 



^ The heating rate could be further corrected by using ~ 
HciB/AnTg^^ when the protons are superconducting. This would 
make core heating even less important at strong fields. 



Ts°° (t) for magnetars are just those with crust heating (see 
Figs.[5]and[7l respectively), and crust heating does not affect 
the onset of core neutron superfluidity because of thermal 
decoupling between the outer crust and core. 

Magn etar activity may be driven by fl eld decay in 
the core (jThompson fc DuncanI 19951, 199(: ) oi_by pro- 
cesses in the crust (iThompson fc DuncanI Il993l . Il995l : 
iGlampedakis et all I2OIII : IPrice et all I2OI2I I. Although the 
former may still be true, it must be accompanied by heat- 
ing in the (outer) crust. On the other hand, fleld evolution 
in the crust easily couples to surface emission. Thus in or- 
der to understand the high surface temperature of magne- 
tars, detailed studies should focus on ma gnetic fleld evolu- 
tion and heating in the crust (see, e.g.. [ Pons et al.l l2009l : 
ICooper fc KaplanI I2OI0I : IPrice et al.1 l2012f ). We note that 
seve ral magnetars may h ave very similar X-ray luminosi- 
ties (jPurant fc van Kerkw iik 20Qg); this could suggest that 
their crustal field strengths are similar and the field decay 
timescale is longer than the age of the oldest of these sources. 

Finally, the (normal versus superfiuid) state of the core 
has important conseq uences for magneti c and rotational evo- 
lution of magnetars (|Glampedakis fc A n dcrsson 2011, ), as 
well as their glitching behavior and possible stellar oscilla- 
tions; studies of these effects may be effective probes of the 
neutron star core. We showed that the core can be treated 
as being in a superfiuid and superconducting state after the 
neutron star is a few hundred years old. 
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